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ABSTRACT 

Context. Applying photometric catalogs to the study of the population of the Galaxy is obscured by the impossibility to map directly 
photometric colors into astrophysical parameters. Most of all-sky catalogs like ASCC or 2MASS are based upon broad-band photo- 
metric systems, and the use of broad photometric bands complicates the determination of the astrophysical parameters for individual 
stars. 

Aims. This paper presents an algorithm for determining stellar astrophysical parameters (effective temperature, gravity and metallic- 
ity) from broad-band photometry even in the presence of interstellar reddening. This method suits the combination of narrow bands 
as well. 

Methods. We applied the method of interval-cluster analysis to finding stellar astrophysical parameters based on the newest Kurucz 
models calibrated with the use of a compiled catalog of stellar parameters. 

Results. Our new method of determining astrophysical parameters allows all possible solutions to be located in the effective 
temperature-gravity-metallicity space for the star and selection of the most probable solution. 
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1. Introduction 

' Estimating the astrophysical parameters of stars is one of the 
, core problems in astrophysics and the key solving a number of 
' astrophysical questions. Indeed, the knowledge of effective tem- 
perature, luminosity class (or gravity at the surface of the star), 
and metallicity of the star allows the distance to the star, the age 
' of the star, and the extinction in the interstellar medium between 
\ the star and the observer to be estimated. 

The number of stars with detailed spectral information is 
very limited compared to the wealth of data in photometry. The 
1 problem of determining astrophysical parameters with photom- 
etry is usually solved with the help of the narrow-band photo- 
metric systems (for instance, the Stromgren photometric system) 
that provide quite good precision for a temperature, gravity, and 
metalhcity. For example, Nordstrom ( 2004) estimates the metal- 
licity from Stromgren photometry with the final average 0.12 
dex dispersion around the spectroscopic values. The problem be- 
comes much more complicated in the case of broad-band pho- 
tometry where all fine features in the photometric bands (tracers 
of the differences in metallicity and especially in gravity) are ef- 
fectively smoothed. The estimation of 7"c//, logg, and [M/H] is 
crucial for the study of the Galaxy on the basis of broad-band 
all-sky surveys like 2MASS or the broad-band deep surveys like 
SDSS/SEGUE or UKIDSS. 

The history of stellar parameter determination from broad- 
band photometry was started by Johnson in the middle of the 50s 
of the last century (see, for example, Johnson |1953i where the 
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author defined extinction-free colors for the UBV photometric 
system). Later, some broad-band photometric systems showed 
the usefulness of extinction-free colors for determining of tem- 
peratures and gravities with some assumptions about the metal- 
licity (for example, the Vilnus photometric system). The use of 
extinction-free colors is discussed in detail in Straizys ( I1977I I. 

A major problem of using extinction-free colors is how to 
determine extinction-free colors themselves from the observed 
ones. Indeed, as we show below, to calculate the extinction-free 
indices and to use them for the estimation of astrophysical pa- 
rameters we have to accept some initial guesses about the esti- 
mated astrophysical parameters of the star a priori. 

In this study we minimize the amount of information to be 
assumed a priori for the estimation of astrophysical parameters 
of the star. We will not make guesses about the astrophysical 
parameters of the star (temperature, gravity or metallicity) based 
on the position of the star, e.g. relative to the galactic plane. 

We do not discuss the dependence of the estimated parame- 
ter on the synthetic model grid adopted in this study and do not 
compare diff'erent synthetic model grids. The estimated parame- 
ters will always depend on the consistency of the synthetic mod- 
els used in the study. However, in the case of Johnson B,V,J,H, 
and Ks photometry, we calibrated the parameters obtained from 
the model fluxes by our method with published parameters from 
different authors. 



2. The synthetic grid of magnitudes and colors 

To start with the determination of the astrophysical parameters 
from photometry we have to provide a synthetic grid fine enough 
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Fig. 1. The synthetic grid of stellar atmospheres used in the 
study. Open circles are original models calculated by Castelli 
(see Castelli & Kurucz'2003 ) and the closed circles are models 
calculated by us with the use of ATLAS9 program to fill the gap 
in the original grid. 

to give us synthetic magnitudes and colors for each point of the 
grid. 

In this section we describe the definition of terms used in 
the paper, the computation of the synthetic grid of colors and 
magnitudes, the influence of extinction and the definition of the 
extinction-free colors. We create a grid of synthetic colors, use 
an assumed extinction law (Cardelli et al. 1 19891 ) to find reddened 
colors and to construct finally a grid of a synthetic extinction- 
free indices. 
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Fig. 2. Comparison of radii used in this study with the radii 
for the main sequence (diamonds), giants (triangles), and super- 
giants (squares) from Schmidt-Kaler ( 1982). The solid curve is 
the theoretical radius-temperature dependence for the main se- 
quence, the dotted curve the theoretical radius-temperature de- 
pendence for giants, and the dashed curve the theoretical radius- 
temperature dependence for supergiants. 



2.1. Definition of syntinetic magnitudes 

The definition of synthetic magnitudes depends on the type of 
detector used in the observation that we are going to approxi- 
mate. If the detected flux corrected for the atmospheric extinc- 
tion is and the stellar flux at the surface of the star is (this 
is the flux delivered by synthetic atmosphere models), they are 
related to each other as 

fA ^ 10-'^'^' (R/dfFi (1) 

where is the interstellar extinction at the wavelength A, R the 
stellar radius, and d the distance to the star, R and d are in the 
same units. In the case of photon-counting devices the magni- 
tude of the star nis in the filter S a is 

f^'Af^SAdA 

ms = -2.5 log + m° . (2) 

l,;ASAdA 

We discuss the zero-points of photometric systems ot° in Sect. [3] 

2.2. Tlieoreticai models used in tine study 

We used the latest Kurucz models of theoretical spectra (Castelli 
& Kurucz 120031 ). because these models can provide wider cov- 
erage in parameter space compared to other libraries. The grid 
used in the study consists of models with temperatures from 
50,000 K to 3,500 K, log^ from 0.0 to 5.0, and metallicities 
from [M/H]= -2.5 to 0.5. We extended this initial grid of mod- 
els to lower gravity and added two points on the gravity axis 
(logg = 2.3 and logg - 2.7) using the program ATLAS9 by 
Castelli & Kurucz (2003). The final grid of original models and 
newly calculated models is presented in Fig[T] These models do 
not have without overshooting, but do have microturbulent ve- 
locity ^ - 2 km/s. The grid was interpolated with cubic splines 



to provide a step in A log T^ff - 0.005 dex, a step A logg = 0.1 
dex, and a step A [M/H] = 0.1 dex. 

To calculate absolute magnitudes from Kurucz models, we 
had to adopt stellar radii for each point in the r^,// - log g grid. 
We took radii from the theoretical track systems of Palla & 
Stahler (1993), D'Antona & Mazzitelli (1994) for pre-main se- 
quence stars, and Geneva tracks (Schaller et al. 119921 ) for post- 
main sequence stars. Figure |2] shows the agreement between the 
radii from the compiled track system (points) and the "stan- 
dard" radii for main sequence stars, giants, and supergiants from 
Table 23 of Landolt-Bornstein (Schmidt-Kaler ,1982 ). The coin- 
cidence is generally quite satisfactory except in the case of hot 
(log Teff > 4.3) supergiants. 

We used the following sets of filters in the study: Johnson 
B and V filters from the UBVBUSER code (Buser [19781 ) and 
2MASS J, H, and Ks filters (Cutri et al. |2006b . 

2.3. Extinction in fhe case of a broad-band pliotometric 
system 

To compute a set of extinction-free indices, we have to compute 
first reddened magnitudes and colors for a wide range of inter- 
stellar extinction. This gives us synthetic color excesses for each 
model in the grid. We adopt the extinction law Aa from Cardelli 
et al.( |1989l ). and take Ay = 3.11 Eb-v, computing a grid of mod- 
els with the step in extinction AAy - 0.05 mag. 

The extinction law used in the study (Cardelli et al. 1 19891 ) 
was constructed with the use of O- and early B-stars. This extinc- 
tion law is a polynomial approximation of the extinction for the 
wavelength /I < Aa/Av >= a(l / A)+b(l / A)/Rv. In fact, this value 
of Ry is a parameter related to OB stars. This also means that 
extinction constructed with the use of an extinction law Aa/Av 
will give different values of the extinction for different spectral 
classes in the same photometric band and with the same dis- 
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Fig. 3. Differences in the extinction law. Solid lines are the differ- 
ences in true and nominal extinction in the Johnson V band for 
stars of different temperatures for different gravities (from log g 
0.0 to 5.0 with 0.5 dex step size). Solar metallicity is assumed. 



tance and ISM between star and observer Let us estimate this 
difference in extinction in the Johnson V band depending on 
the spectral class and luminosity class of the star. We calculate 
the extinction as the difference between theoretical magnitudes 
for the Johnson V filter and the stellar atmosphere model with 
logT,ff,logg,[M/H] as A"^"''' = my(A7"') - myCA^"' = 0) 
with the use of A^(Av) from Cardelli et al. ( 119891 ): 



^model ^ _2_5 log 



A fx SxdA 



(3) 



The output Ay"''''' differs from the input A'^''"' and this difference 
depends on the spectral class, gravity, and metallicity of the star 
This difference is a physical phenomenon. Indeed, an 05 star and 
an M5 star being placed at the same distance from the observer 
and with the same optical depth of the ISM will be differently 
obscured by the ISM. This difference reaches up to 5% of the 
nominal value of an extinction for the later spectral types. 



2.4. Extinction-free color indices 

Extinction-free color indices (or Q-values) are the final target 
for our computations of synthetic colors. With the use of these 
indices, we are able to reduce the number of models in our grid, 
because we can use a unique grid of synthetic models for any 
star independent of the individual extinction for this star. 

The construction of extinction-free color indices has some 
problems that we consider below. The Q-values are by definition 



Qm\n 



- (nil - m2)o - 

- (nil - m) - '■ 



(m2 - m3)o 
■ (m2 - m-i). 



(4) 



However, Em,,n2/Em2my is not constant, but depends on the spec- 
tral class, gravity and metallicity of star (see, for example, 
Grebel & Roberts 1995 ). The excess ratio depends as well on 
the extinction itself. By definition, £■„,,, „j = A,„, - A,,,^, so to 
remove extinction from the excess ratio, we have to replace 
A„,, with A,„, or vice versa. However, the dependence of the 



extinction on the optical depth is not linear (see, for instance. 
Chapter 6.5 from Straizys [19771 1. As a result, we have a depen- 
dence A,„, - ki A„„ -I- ^2^m2 + • ■ -1 and to construct an excess 
ratio, we have to neglect orders higher than the first in this de- 
pendence. For the precise determination of astrophysical param- 
eters, we have to check the possibility of neglecting these higher 
order coefficients. Figure |4] shows the remaining dependence of 
the Q-values on extinction, which can indeed be neglected com- 
pared to the Q-values themselves in Fig.|5] 



3. Calibration of theoretical models 

As we see from Eq. (|2]i, we have to bring our theoretical models 
into accordance with the zeropoint of the photometric system 
(mO). 

To provide this calibration, we compiled a catalog of stars 
with temperature, gravity, and metallicity known from spec- 
troscopy. The catalog is a combination from the catalogs of 
Borkova & Marsakov (120051 ) and Cayrel de Strobel et al. ( 1997ll 
with the addition of the catalogs of Soubiran el al. (1998), 
Erspamer & North ( 2003) and Takeda et al. (2005). Our final list 
contains 3092 stars, 2096 of them were identified with the ASCC 
catalog (Khai-chenko'2001V The latter contains 2,501,304 stars 
with Johnson B, V photometry recalculated from the Tycho Bj, 
Vt and parallaxes and proper motions (Hipparcos and Tycho-1). 
We crossmatched the ASCC catalog with the 2MASS catalog, so 
that the stars in our compilation (2096 stars which we will call 
"calibration catalog" further on) have Johnson B, V and 2MASS 
J, H, Ks magnitudes. 

A subsample of 1287 stars was selected with ctb < 0.1 mag, 
cry < 0.1 mag, cry < 0.3 mag, cth < 0.3 mag, ctk, < 0.3 mag, 
parallaxes with <T„/n better than 10% and distances within 100 
pc. We assume that the stars from this subsample are extinction- 
free and use the subsample to calibrate our theoretical models. 

To calibrate the theoretical sequences, we selected all stars 
from the subsample with log^ = 4.0 + 0.1, [M/H] - 0.0 + 0.2, 
and photometric errors in all bands better than 0.05 mag and 
absolute magnitudes in the B band Mb < 6.0 mag. 

On the basis of this calibration subsample of 58 selected 
stars, we calibrated the synthetic models 



(B - Vf = (B- Vf + (0.6001 + 0.0445) 

(V - jf = (V- Jf + (2.6764 + 0.0896) 

(J-Hf^ (7-//)^ + (1.0382 + 0.0361) 

(H-Ksf ^ (H-Ksf + (1.1141 + 0.0290). 



(5) 



Although, the parameters from literature cover only a small frac- 
tion of the full parameter space, we applied the corrections to 
the full sythetic grid. Comparison of our calibrated colors with 
the colors for Vega shows good agreement. We take the atmo- 
spheric parameters for Vega from Smith & Dworetsky ( 119931 ) 
(T,ff = 94507i:, logg = 4.00, [M/H] =-0.55), which from our 
calibrated model grid yields (B - V) = 0.005 mag, (V - J) = 
0.014 mag, (J - H) = -0.060 mag, (H - K,) = 0.052 mag. 
The 2MASS JHKs photometric system differs slightly from the 
Johnson JHK (see Carpenter i2001i and Koen et al.|2007j). 

4. The interval-cluster method of the astrophysical 
parameters estimation 

Parameters like effective temperature, gravity, metallicity, ex- 
tinction, and distance on the basis of broad-band photometry 
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Fig. 4. Dependence of the Q values on Ay for stars with solar 
metallicities, \ogg = 4.5, and Te/f = 3500 K (filled circles), 
Teff = 7000 K (open circles), and Te/f = 10000 K (crosses). 



must be determined directly from photometry without any a pri- 
ori assumption about physical properties of stars that can restrict 
the range of the search. Suppose that we have a parameter space 
{XQ}={Teff,logg, [M/H], Ay, r), where T,,// is the effective tem- 
perature, \ogg the logarithm of gravity, [M/H] the metallicity. 
Ay the extinction in the Johnson V band (Ay is used as the pa- 
rameter for extinction hereafter), and r the distance to the star. 
From observations (or from a synthetic grid in the case of simu- 
lations), we have a magnitude space {Y^)]-{mi, . . . , m„] that can 
be transformed to the color space { Y ]-{ m, - nij, i + j]. Via 
this transformation, we eliminate the distance r from the param- 
eter space {Xq] and work in the parameter space {X]={Teff, log g, 
[M/H], AW. 

Let us suppose that there are two different subsets in the pa- 
rameter space { X^ } and { X^]. These two classes can be popu- 
lated by stars within any given temperature range, gravity range, 
and metallicity range. Let us as well have two subsets in the 
color space {F'} and {Y-] containing the synthetic colors of the 
elements in {X' } and [X^]. The task of classification is to attribute 
the star with some colors to one of the classes in the parameter 
space (X}. This is possible if the two subsets [Y^] and [Y^] are 
separated in the color space, but this requirement is not always 
guaranteed in the case of broad-band photometry. Stars with dif- 
ferent temperatures, gravities, metallicities, and extinction can 
occupy the same area in the color space. 

Figure |6] illustrates the problem: stars within two different 
metallicities and gravities ([M/H]=-0.5, log^ € [4.5,5.0] and 
[M/H]=0.0, \ogg e [1.0,2.0]) occupy the same space in the 2- 
color diagram even in the absence of extinction. It is obvious that 
for this range of colors there are at least two different solutions 
in parameter space for the star inside the filled area. 

In this case, there are at least two possible solutions for the 
parameters of a star: { X^ : [M/H] = -0.5,logg 6 [4.5,5.0] } 
(dwarf) and {X^ : [M/H] = O.0,logg e [1.0,2.0]) (supergiant). 
The idea of the method we propose here is the location of all 
probable solutions {X'] and an estimation of the probability that 
the star can belong to the solution {X']. 
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Fig. 5. The Q diagrams. Panel (a) shows the Qbvj-Qvjh diagram 
for 3 different metallicities (the thick line is for [M/H]=0.5, the 
intermediate line for [M/H] =0.0, and the thin line for [M/H]=- 
0.5), and \ogg - 4.5. Effective temperatures for [M/H]=0.5 are 
shown, effective temperatures for other metallicities are, respec- 
tively, 50000 K for the point with maximum Qvjh, 4300 K for 
the point with minimum Qvjh, and 3500 K for the end point 
of the sequence. Panel (b) is the Qbvj-Qvjh diagram for so- 
lar metallicity and \ogg 0.0 (dotted-dashed line), 3.0 (dashed 
line), and 4.5 (solid line). Effective temperatures are shown for 
log g -0.0 (in italics, 7000 K for the point with maximum Qv jh, 
3900 K for the point with maximum Qbvj and 3500 K for the 
end point of the sequence), \ogg = 3.0 (30000 K for the point 
with maximum Qv jh, 4000 K for the point with minimum Qv jh, 
and 3500 K for the end point of the sequence), log g = 4.5 (in 
bold, 50000 K for the point with maximum Qv jh, 4300 K for the 
point with minimum Qv jh, and 3500 K for the end point of the 
sequence). Panels (c) and (d) are the same but for Qbvj-QjhKs 
diagram. 

To reduce the number of parameters in the parameter space, 
we transform the colors of the star into the extinction-free color 
indices described above. The method is illustrated in the follow- 
ing 5 bands: Johnson B and V magnitudes (B,V filters are from 
UBVBUSER code, Buser 19781 and 2MASS J, H, and Ks mag- 
nitudes (Cutri el al. 2006). Four colors {B - V, V - J, J - H, 
H-Ks) and three extinction-free color indices (Qbvj, Qv jh, and 
Qjhks ) are used. 

Our method consists of three steps: 

1. interval analysis: select all theoretical models in our grid 
which satisfy the extinction-free color indices of the star; 

2. cluster analysis: group all possible solutions in parameter 
space and find a set of solutions; 

3. selection of the best solution: consider all possible solutions 
to select the only one as the most probable solution for the 
star in process. 

The final solution for the star is a set of parameters (mean 
values and variances): {< Tgff >', < logg >', < [M/H] >', 
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4.5 and 5.0, dashed lines for solar metallicity and log g between 
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< Av >• 



/el, «, where n is the num- 



ber of possible solutions for the star. One of these solutions will 
be selected as "the most probable" due to criteria we describe 
below. 



4. 1 . The first part of the method: interval analysis 

In this step we locate possible solutions in parameter space 
(log r^yy, log^, [M/H]). The table of the theoretical models is 
constructed with steps AlogTe// - 0.005 dex, Alogg = 0.1 
dex, and A [M/H] = 0.1 dex and has 249922 entries in the 
parameter space. For each entry there is one combination of 
theoretical colors {{B - V^AV - Jf,{J - HY AH - Ks f ) 
and extinction-free color indices {{Qbvj)tXQvjh)tXQjhk)t), as 
well as the set of supplementary parameters: ratios of extinction 
((Ab/AvV, (Aj/Avf, (Ah/AvV, (AKjAvf), ratios of redden- 
ing {{Ebv/EvjV, {Evj/EjhV, {EjhIEhksY). The "observed" 
values for the star are magnitudes (B^ ,K^) and colors 

((B - Vf, (V - Jf, (J - Hf, (H - K)l). 



The difference between synthetic Ql 



(B - vy 



(Ebv/EvjV (V-Jf and the "observed" index Qlyj = (B-Vf - 
{EBvlEvjf {V-Jf in this ideal case will be zero. However, we 
know nothing about parameters of the star except our assump- 
tion that the stellar parameters are inside our synthetic grid. 

Let us try each of the grid points for the input colors of the 
star. We take the difference 



AGfivy = \Qlyj -{B-Vf + {V- 



(7) 



ivy 



in each grid point. We do the same for AQvjh and AQjhks ■ All 
points that satisfy the following equations 



^Qbvj < cr Qgyj 
^QvjH < cr Qyj„ 

are selected as possible solutions. 

The definition of (tq^^, comes from the definition of Q: 



(8) 



Qbvj 



(B-vy 



+ (■ 



■'BV 



(9) 



The problem are the quantities of {Ebv I EyjY and cr^ 



which are unknown a priori. In this step we take the maximum 
value of {EbvIEvjY in our theoretical grid (averaged over all 
points in the grid < (Ebv/EvjV >- 0.426466 + 0.156561, max- 
imum value max (EbvIEvjY - 0.4516) for all points, and dou- 



ble the value of (Ebv/Evj ) 
Finally, 



cr. 



Qbvj 



= cr, 



(B-V)s 



■ 0.9032 (rly_j^^. 



as a rough estimation of 



(10) 



The same is done for the estimation of o-q^j^ and o-Qj^^. . 

After the initial selection of grid points in accordance with 
Eq. (O we return from the extinction-free space to the color 
space and remove all "artifact" entries (those with negative ex- 
tinction and colors out of the bounding box for the input star). 

We assume that the ratio Ab/Av is constant for a given tem- 
perature, gravity, and metallicity and that is only depends on 
these parameters and not on Ay itself. 

For each of the selected grid points, we can find reddening 
for each color and extinction for each band. 



E/ ^c' -J 



(11) 



where is the observed color for the star and is the synthetic 
color at the selected grid point, c = (B - V), (V - J), (J - H), (H - 
Ks). Also, 



4.1 .1 . Selection of grid points 

The selection of grid points satisfying the set of "observed" col- 
ors is carried out in the space of extinction-free color indices. 
The value {QbvjV does not depend on extinction. If we know 
stellar parameters (log ff, log g, [M/H]), we can find a stellar 
excess ratio Ebv/Ev/ ■ Then we can construct an "observed" 
extinction-free color index Qgyj and compare it with the the- 
oretical ones Q^yj to select grid points satisfying to the input 
index Q: 



E,. 



(12) 



where (AyjEcY is the synthetic ratio of extinction in the 
Johnson V band to color excess in c color at each grid point. 

We use the inverse of the variances of the input colors as 
weights to find the weighted average for the extinction. From 
the previous equations and the definition of Ec^ the weights are 



{EJAy)\ 



(13) 



Qlvj ^{B-Vf-^\v-jf^{B-Vf- — 



Ebv^ ,,5 



(V-7)\(6) 



To remove negative Ay, only grid points with Ay > -3 cr^,, were 
selected, where cr^,, is a weighted dispersion of Ay. 
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Table 1. The calibration subsample 



Catalog 


Number of 
stars 


Interval in 


Tefi 


,K 


log g, dex 


[M/H],dex 


Borkova & Marsakov (|2005) 


876 


4520- 


-6790 


2.94 — 5.00 


-3.85 — 0.55 


Cayrel de Strobel et al. (1997 ) 


3247 


2291 ~ 


36000 


0.09 — 6.00 


-5.6 — 2.9 


Soubiran et al. (1998.) 


211 


3922- 


-6276 


0.58 — 4.69 


-2.91—0.34 


Erspamer & North ll2003ll 


140 


5902 — 


10375 


3.04 — 4.36 


-1.08—1.65 


Takeda et al. d2005b 


160 


5009.3 - 


- 6967.7 


3.185—4.865 


-1.291 —0.457 



For each point in the grid we have theoretical magnitudes, 
colors, theoretical excess ratios, and an estimation of the extinc- 
tion Ay- Now we calculate for each point a synthetic reddened 
color. 



cl ^c'' + Av, 

Ay 



(14) 



and remove points that are out of the 3 cr range from the input 
color of the star. Finally, we select points with 



\cl -c^\ < 3o-,, 



(15) 



for each color c and have N grid points as possible solutions for 
the star. 



4.2. The second part of the method: cluster analysis 

We selected N points that satisfy our requirement (Eq. |8]l, that 
are within the allowed extinction range {Ay > -3 cr^,,), and that 
are inside the error box of the stellar colors. All these points are 
possible solutions for the input star must be grouped properly to 
form a set of possible solutions in the next step. 

The distribution of the grid points is uniform in the parameter 
space (log Tgff, \ogg, [M/H]). We selected all grid points that are 
inside the error box of stellar colors, so are now able to simulate 
a more realistic physical distribution of grid points inside the 
error box assuming a normal distribution of errors of the colors 
of the star This step will be made by giving new weights to each 
selected gridpoint. 

The weighted distribution of grid points inside the error box 
will be projected onto the axes in the parameter space to subdi- 
vide the subspace with selected grid points into the set of possi- 
ble solutions and to assign a probability to each solution. 

4.2.1. Weighting of points 

Suppose that the "observed" star with colors {B - V)^,(V - 
J)^ ,{J - H)^ ,{H - Ks)^ has a normal distribution of errors in 
each color (A^(0, cr^)) and all colors are independent variables. 
From Eq. (fT4l i we can find the distance of each grid point from 
the center of the error box (from the observed color of star). 



Ac : 



(16) 



where c is one of the colors, c = (B-V), (V-J), (J-H), (H-Ks). 

The probability that a point occupies the 4-dimensional color 
space from (B-V) to (B-V)+d(B-V), (V-J) to (V-J)+d(V-J), 
(J-H) + d(J - H), (H - Ks) + d(H - Ks) is 



dN 



n 



/(Ac)flfc, 



C=(B-V), (V-J), (J-H), {H-Ks ) 



where 



f(MB -V))^e 



(17) 



(18) 




3.56 3.60 3.64 

og Teff 



3.68 



Fig. 7. The projection of the weighted distribution function of 
selected grid points on the log Tefj axis, showing the subdivision 
made with an assumed 10% significance level. 




3.56 3,60 3.64 3,68 
log Teff 



.56 3,60 3.64 3,68 
log Teff 



Fig. 8. Points selected for the test input point in the parame- 
ter space. Big open circles are for solution 1 (see Table |2]i, big 
closed circles for solution 2, crosses for solution 3, and small 
circles are for solution 5. 



The same distribution for other colors is assumed. 

The real distribution of grid points is N( A(B - V), A( V - 
J), A(J - H), A(H - K)). This is a number of points with A(B - 
V), A(V - J), A( J - H), A(H - K) between A(B - V) and A(B - 
V) + d(B - V), A(y - J) and A(V -J) + d(V - J), and so on; i.e. 
Np = N(A(B - V), A(V - J), A(J - H), A(H - K)) d(B - V) d(V- 
J) d( J - H) d(H - K). 

The one-dimensional weight for the color (B-V), for exam- 
ple, is wa(b_v) = f(A(B - V)) d(B - V) and for a point in the 
4-dimensional color space 



n vvAc 

c={B-V),{V-JUJ-HUH-Ks) 



(19) 
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Table 2. Set of solutions for a test point of T^jj - 4000 K, log g - 4.45 and [M/H]=0.0 and with an assumed error in color of 0.05 
mag. 



Number 
of points 


Total 
weight 


J-' mm 

<■//' 

K 


eff ' 

K 


K 


logg""", 


logg""', 


logg"' + 0-tog^, 


[M///]™" 


[MIHT"' 


[MIHY ± cr^Mim 


1713 


0.7204 


3664 


4305 




2.2 


5.0 


3.8 ± 0.4 


-2.5 


0.5 


-0.9 + 0.8 


951 


0.0744 


4305 


4721 


4515-'" 


0.3 


3.8 


2.2 ± 0.5 


-2.5 


-1.9 


-2.2 ± 0.04 


2736 


0.1346 


4305 


4721 


4561-'*' 


0.9 


5.0 


3.5 ± 0.7 


-1.9 


-0.3 


-1.2 + 0.2 


1 


0.0006 


4305 


4305 


4305 


5.0 


5.0 


5.0 


0.5 


0.5 


0.5 


1728 


0.0700 


4721 


4943 


4754-'' 


0.3 


5.0 


2.6 ± 1.3 


-2.5 


-0.3 


-1.4 + 0.4 
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4 •I'liliF'- 
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en 
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4 
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c 
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5 

4 
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1 


5 

4 

° 2 
1 
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(10) 



+++- 



+-++++ 



(2c) 



3$S 




(3a) 



t . > t M f t ■ H ■ ■ ^4 + 

: * : • -1 * 4- ^ • 

t 



(4a) 



r 



t I I M 



(lb) 



(2b) 



I 



(3b) 



— I ^ * 

t M h ' 
< I I 



(4b) 



3-6 3.8 4.0 4.2 4.4 4.6 3.6 
log T,f, 



3.8 4.0 4.2 4.4 4.6 
log Tan 



Fig. 9. The offset (bias) in the case of the selection of the solu- 
tion with the largest weight in the (log T^jy, logg) plane. Open 
circles are the cases of a single solution, lines show the distance 
between the selected solution and the input grid point. Panel 
(la) is for metallicity [MjH} - -2.0 dex and an assumed rms- 
error cTph - 0.01 mag in all colors, panel (la) for metallicity 
{MjH} - -2.0 dex and an assumed rms-error cTph - 0.05 mag 
in all colors, panel (2a) — [MjH} = -1.0 dex and CTph = 0.01 
mag, panel (2b) — {MjH} - -1.0 dex and cr^/, = 0.05 mag, 
panel (3a) — {MfH} - 0.0 dex and cr^/, = 0.01 mag, panel 
(3b) — [MIH] = 0.0 dex and cr,,/, = 0.05 mag, panel (4a) 
— [MIH] = 0.3 dex and a-ph = 0.01 mag, panel (4b) — 
[MIH] = 0.3 dex and cTph = 0.05 mag. 



The optimal bin in the case of the one-dimensional distribution 
\s h - 3.5 crn '^^ and, in the case of d dimensions, the optimal 
binis/z - 0(«-'^'''+^') (see Scott etal. 2004). After the weighting 
of the points we normalise the weight of all selected points to 1 , 
so that 2, w' = 1. 

4.2.2. The final set of solutions 

The distribution of the weighted grid points was projected on 
each axis in the parameter space (log T^fj, log g, [M/H]) to group 
the grid points into clusters in the parameter space. These clus- 
ters are the final set of possible solutions for the input star. 



Fig. 10. Errors of the selected solution on the (log re//,log g) 
plane. Panel (la) is for metallicity [MjH} - -2.0 dex and an 
assumed error a-ph - 0.01 mag in all colors, panel (la) for metal- 
licity [M///] = -2.0 dex and an assumed error cr^/, = 0.05 mag 
in all colors, panel (2a) — {MjH} - -1.0 dex and a-pi, - 0.01 
mag, panel (2b) — [MjH} - -1.0 dex and cTpk - 0.05 mag, 
panel (3a) — [MjH} - 0.0 dex and a-pi, - 0.01 mag, panel 
(3b) — [MIH] = 0.0 dex and a-ph = 0.05 mag, panel (4a) 
— [MIH] = 0.3 dex and a-ph = 0.01 mag, panel (4b) — 
[MIH] = 0.3 dex and aph = 0.05 mag. 

Let us suppose that the grid points are distributed between a 
minimum log T"™j. and a maximum log 7""^"- We used the orig- 
inal bin of the grid (AlogTe// = 0.005 dex) to construct the 
projection of the distribution of N selected grid points on the 
temperature axis: 

N 

./logr,,,, = Yj "^j' ■ ^^// = ^^//' ^^^^ 

where i runs through 1 to M, and M is the number of bins 
of the grid between log r™ and log T'^JJ, M = (log T'^JJ - 

log r™p/A log T,ff. Note that log J,//, log g and [M/H] are dis- 
crete values of the grid. 

Let us suppose that we have found L local maxima in 
log Tgff space. We estimate the significance of each local max- 
imum in accordance with the maximum likelihood criterion for 
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Fig. 11. The bias in the case of the selection of the solution with 
the largest weight in the (log Teff,[M/H]) plane. Open circles are 
the cases of a single solution, lines show the distance between 
the selected solution and the input grid point. Panel (la) is for 
logg =1.0 dex and an assumed error o-pi, - 0.01 mag in all 
colors, panel (la) for log g - 1.0 dex and an assumed error o-ph - 
0.05 mag in all colors, panel (2a) — log^ = 2.0 dex and cTph - 
0.01 mag, panel (2b) — logg - 2.0 dex and o-pi, = 0.05 mag, 
panel (3a) — log g - 3.0 dex and o-pi, = 0.01 mag, panel (3b) — 
logg = 3.0 dex and a-pi, - 0.05 mag, panel (4a) — logg = 4.5 
dex and o-pi, = 0.01 mag, panel (4b) — logg = 4.5 dex and 
CTpi, = 0.05 mag. 

the significance of the maximum (see Materne'1979), i.e. for the 
k-th local maximum the criteria for significance is 

Xl^-2ln^, (21) 
where 

^ = n^i4g7^„' (22) 

and is the likelihood in the absence of the k-th local maxi- 
mum. 

We use the 10% significance level (x^ > 2.7). If;^'^ < 2.7 the 
maximum is neglected and joined with the neighboring, most 
significant maximum. This procedure decreases the number of 
clusters in log 7^// space to N < L. In each of these clusters, 
we are looking for a clustering on the logg axis, and for each 
cluster on the logg axis, we are looking for the clustering on 
the metallicity axis with the same procedure. For example, there 
are Mog r^ff clusters on the temperature axis, the i-th cluster on 
the temperature axis has N^^^^ clusters on the gravity axis and 

the j-th cluster on the gravity axis has N'^lj/[j-^ clusters on the 
metallicity axis. 

The final number of clusters (and this is the final set of solu- 
tions) is 

J—!-'-' log,? 

N.,.srer.- J] Z^/H]" (23) 

'■=1 ;=i 



Fig. 12. Errors in the case of the selection of the solution with 
the biggest weight on the (log 7'f,//,[M/H]) plane. Panel (la) is 
for logg =1.0 dex and an assumed error o-pi, - 0.01 mag in all 
colors, panel (la) is for logg =1.0 dex and an assumed error 
CTph - 0.05 mag in all colors, panel (2a) — logg = 2.0 dex and 
CTph = 0.01 mag, panel (2b) — logg = 2.0 dex and o-pi, = 0.05 
mag, panel (3a) — logg = 3.0 dex and a pi, - 0.01 mag, panel 
(3b) — logg = 3.0 dex and cr^/, = 0.05 mag, panel (4a) — 
logg = 4.5 dex and o-pi, = 0.01 mag, panel (4b) — logg = 4.5 
dex and cr^/, = 0.05 mag. 

We have distributed N grid points selected over Nci„siers clusters, 
the sum of weights of points inside the cluster gives us the weight 
of each cluster: 

W; = 2wi,w,:/eA?,-, (24) 

where A^^ is the j-th cluster. 

To find average values of temperature, gravity, and metal- 
licity for each cluster, we repeated the procedure of weighting 
described in Sect. 14.2.11 this time inside each cluster, and com- 
pute the weighted average solutions for logTe//, logg, [M/H], 
Ay, and distance r with the dispersion of these values inside the 
cluster. 

4.3. The third part of the method: selection of the solution 

We forme a set covering all possible solutions for a star The final 
and most complicated step is the selection of one solution from 
the set. 

As we see from Fig. |6] there is no possibility of selecting a 
single "true" solution in general. Nevertheless we can select "the 
statistically most probable" solution by the appropriate criteria. 
We would like to underline that, by these criteria, we do not 
use guesses about the astrophysical parameters themselves, so 
the selection at this point does not introduce limitations on the 
parameters. 

There are a number of methods for this selection: we can take 
the solution with the largest weight or with minimal distance 
from the input "observed" color in color space or take one with 
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Fig. 13. The best regions for estimating astrophysical parameters 
on the (logTeff,logg) plane. Biases and errors of the solution 
(see text for the explanation) are in the range: A log T^ff < 0.05, 
o'iogT,,, < 0.05, Alogg < 0.3, cTiogg < 0.3. Panel (la) is for 
metallicity [M/H] = -2.0 dex and an assumed error cTph - 0.01 
mag in all colors, panel (la) for metallicity [M/H] = -2.0 dex 
and an assumed error o-p/, = 0.05 mag in all colors, panel (2a) — 
[M/H] = -1.0 dex and cr^,, = 0.01 mag, panel (2b) — [M/H] = 
-1.0 dex and o-pi, = 0.05 mag, panel (3a) — [M/H] = 0.0 dex 
and (Tpf, - 0.01 mag, panel (3b) — [M/H] - 0.0 dex and cTpf, - 
0.05 mag, panel (4a) — [M/H] = 0.3 dex and o-pi, = 0.01 mag, 
panel (4b) — [M/H] = 0.3 dex and o-ph = 0.05 mag. 




Fig. 14. The best regions for estimating astrophysical parameters 
on the (log Tgff, [M/H]) plane. Biases and errors of the solution 
(see text for the explanation) are in the range: A log Te/j- < 0.05, 
o-iogT,ff < 0.05, A[M/H] < 0.3, o-[m/h] < 0.3.Panel (la) is for 
\ogg =1.0 dex and an assumed error cTph = 0.01 mag in all 
colors, panel (la) for log g = l.Q dex and an assumed error o-ph = 
0.05 mag in all colors, panel (2a) — logg = 2.0 dex and o-ph = 
0.01 mag, panel (2b) — logg = 2.0 dex and o-ph - 0.05 mag, 
panel (3a) — log^ = 3.0 dex and o-ph = 0.01 mag, panel (3b) — 
\ogg = 3.0 dex and o-ph = 0.05 mag, panel (4a) — logg - 4.5 
dex and o-ph - 0.01 mag, panel (4b) — iogg = 4.5 dex and 
o-ph = 0.05 mag. 



an unbiased distribution of grid points around the center of the 
error box in color space. 

Let us show an example of the determination of parameters 
to point out the major troubles we have to face. We used the 
point from the synthetic grid with Tg/f - 4000 K, \ogg =4.45, 
and [M/H] =0.0, photometric errors equal to 0.05 mag in all col- 
ors were taken. The first step of the method returned 6000 points 
distributed between 3664 K and 4943 K in temperature, between 
0.3 and 5.0 dex in log g and between -2.5 and 0.5 dex in metallic- 
ity. All points were weighted and grouped in 5 possible solutions 
(see Table |2]i. The second part of the method (the cluster anal- 
ysis) is illustrated by Fig.|7] The first maximum on the log Teff 
axis is rejected by the significance level assumed. Finally Fig. [8] 
presents the clusters in the projections onto the log T^ff - log g 
and the log Tg/j- - [M/H] planes. 

We define the following two quantities that describe the error 
of the estimated parameter: 

1. The offset (bias) of the selected solution from the true pa- 
rameters AlogTeff, Alogg, A[M/H]. These values are the 
differences between the estimated parameter (the center of 
the selected cluster) and the true parameter known only in 
simulations. 

2. The dispersion of the estimated parameter o-iogT^ff, o-\ogg, 

O" [M/H] ■ 



From Fig. [8] we see that it is possible to subdivide the first 
cluster (first solution in Table |2]i. This possibility was ruled out 
by the significance level we assumed (see Fig.|7] the first peak on 
the distribution function). If we decrease the significance level, 
we are increasing the number of clusters in the set of solutions 
but decrease the weight for each cluster Selecting a single clus- 
ter we face the problem of choosing between a number of clus- 
ters: the one that includes the true parameters of a star we call 
"true" cluster (solution) for the star, the remaining are "false" 
solutions. In practice we can distinguish between "true" and 
"false" solutions only if we knew the parameters of the star a 
priori; however, we can determine the region in the parameter 
space where only a single, hence "true", solution is available or 
where the "true" solution can be selected from a set of solutions 
by some criteria. We note that the definition of this region in 
parameter space depends on the desired precision of the deter- 
mination of parameters. We show below how to define such a 
"most preferable" region. 

4.3.1. The most preferable regions in parameter space 

We have determined regions in the parameter space where we 
can best estimate astrophysical parameters; in other words, if the 
solution falls into this region, we can be sure that this estimation 
is a correct, unbiased solution for the star. 
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Fig. 15. The test of the method on the cahbration sample estimating temperature. In the left panel there are all stars with cr\ogT,ff < 
0.05 for an estimated temperature and errors of photometry better than 0.05 mag in all colors. The middle panel shows a subsample 
of the left panel falling into the most suitable region on the (log T^ff, \ogg) plane as described in the text and the right panel gives a 
subsample of the left panel v^ith Ay[mag] I d[kpc\ e [0.1,2.0]. The dashed lines confine the 3-cr region; i.e, they are the median line 
and log r/ff = log rf^^ + 0.15, log T^'l^ = log rj^ - 0.15 lines. 
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Fig. 16. The test of the method on the cahbration sample for estimating gravity. In the left panel there are all stars with criog^ < 0.5 
for the estimated gravity and photometric errors better than 0.05 mag in all colors. The middle panel shows a subsample of the left 

panel falling into the most suitable region on the (log Tejs^ log^) plane as described in the text, and in the right panel, we plot a 
subsample of the left panel with Aylmagj/dlkpc] e [0.1,2.0]. The dashed lines are the median line and log^*^ = log^'''' + 1.5, 
log^*^'=log^''''-1.5Unes. 



Table 3. Selection of the solution: comparison of different cases 



Parameter 


Number of stars 


average and variance 




criterion 


criterion 


criterion 


criterion 


criterion 


criterion 




1 


1 &2 


1 &3 


1 


1 &2 


1 &3 


log r;;^- log r;" 

logg'P-log^;"'' 


588 


31 


330 


0.001 ± 0.059 


0.054 ±0.176 


-0.010 ± 0.025 


133 


38 


40 


-0.166 ±0.924 


-0.336 ± 1.047 


0.191 ±0.778 


[M/HY" - [M/H]"'' 


527 


199 


305 


0.011 ±0.457 


-0.029 ± 0.459 


0.017 ± 0.407 



To find the most preferable regions, we performed simula- 
tions using the synthetic grid points as input "observed" col- 
ors. We interpolated between the grid points to create from 
a grid point with {logTgff, logg, [M/H]) new points with { 
\ogTeff+iA\ogTeff/2, logg+jMogg/2, [M/H]+kA[M/H]/2 
}, i,j,k = [0, 1], and A log 1^//, A logg, and A[M/H] are inter- 



vals between points of the initial grid in temperature, gravity, and 
metallicity. For each point of this enhanced grid, we calculated 
synthetic colors. The resulting test grid includes points with the 
largest scatter from the points of the initial grid in the parameter 
space. 
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Fig. 17. The test of the method on the calibration sample for estimating metallicity. In the left panel there are all stars with cryMiH] < 
0.3 for the estimated metallicity and photometric errors better than 0.05 mag in all colors. The middle panel shows a subsample of the 
left panel falling into the most suitable region on the (log T^fj, [MjH]) plane as described in the text, and in the right panel, we plot 
a subsample of the left panel with Av[mag] /d[kpc] £[0.1, 2.0] . The dashed lines are the median line and [M/H] - [M/HY'' + 0.9, 
[M/HyP = [M/HV'' - 0.9 lines. 



We tested the method on two sets of errors in colors; in the 
first case, all colors had an error 0.01 mag and in the second 
all colors had an error 0.05 mag. Figure |9] shows the solution 
selected by the maximum weight with offsets (biases) caused 
by the selection of the "false" cluster or the shift inside the se- 
lected "true" cluster, and Fig.[TT|is the same for the temperature- 
metallicity plane. Figure [10] plots the rms-errors of the selected 
solution for the temperature - gravity plane and Fig. [12]- for the 
temperature-metallicity plane. 

We define a region in the parameter space as suitable for the 
unbiased determination of astrophysical parameters, if we can 
select a solution for the input point inside the region, and this 
final solution has biases for temperature AlogTe/f < 0.05, for 
gravity Alog^ < 0.3, rms-errors criogT.jj < 0.05, criogg < 0.3 in 
the temperature-gravity plane and A log T^ff < 0.05, A[M/H] < 
0.3, rms-errors o-iggj-,/! < 0.05, o-[m/h] < 0.3 in the temperature- 
metallicity plane. This is done in Figs. [13] and [14] According to 
these figures we define the following regions in the log Te/f - 
log g plane 

log r,//< 3.75, logg> 4.0; 
logT.ff e [3.395, 3.7], logg < 1.0; 
logTeff e [3.8,4.05], log g > 4.5; 
logr,//G[3.9,4.2],log^e[2,3]; 
log Teff > 4.4, \ogg > 4.0 

as the most suitable for the determination of log T^jf and logg, 
and on the log T^ff - [M/H] plane 

logr,// 6 [3.7,3.8] 

as the most suitable in the search for metallicity. We note that 
the definition of these regions depends on the assumed errors of 
photometry, i.e. observations; and in the case of real stars with 
different observational errors in different bands, the definition of 
the most suitable regions can vary. 

4.3.2. Setting bounds on extinction 

Estimating the most preferable regions for our method made in 
the previous section is based on the assumption that photometric 



errors are the same in all bands. In the case of real observations, 
we would have to recalculate the estimation of errors and biases 
at any point on the synthetic grid for each combination of er- 
rors of colors for the input stars. This is possible but consumes 
enormous computational resources and slows the estimation. 

For the actual application of the method to observations, we 
propose the following procedure: we estimate the extinction per 
photometric distance. We assume a maximum extinction in the 
galactic plane to be 2.0 mag/kpc in the Johnson V band and a 
minimum extinction in the direction of the galactic poles as 0. 1 
mag/kpc (from Schlegel et al. |1998l ). 

Finally there are three criteria for selecting a single solution 
for the star: 

1 . selection of the solution with the largest weight, 

2. selection based on the most preferable region in the parame- 
ter space, 

3. selection based on setting a bound on the extinction per pho- 
tometric distance. 

5. The method used in practice and problems of the 
method 

We used the 2096 stars of our compiled catalog (see Sect. [3] for 
details) with Johnson B, V, and 2MASS /, H, Ks photometry to 
test the method in practice. 

Table|3]shows the statistics for the different combinations of 
the criteria described in the previous section: criterion 1, crite- 
rion 1 combined with criterion 2 for selected stars (we select a 
single solution with the largest weight for each star and select 
stars within the most preferable regions), and criterion 1 com- 
bined with criterion 3 for selected stars. The last 3 columns in 
the table show the means and variances of the differences be- 
tween the "true" spectroscopic parameter and our estimation (for 
example, < log r'j!^ - log T'fl^ > ±(^iogri^-\ogT'jy The results 
for the calibration subset, based on the different selection criteria 
are shown in Figs. [15] [H] [17] (left panel for criterion 1, middle 
panel for criterion 1 & 2, right panel for criterion 1 & 3). 

The results for all stars in the compiled catalogue are shown 
in Figures [18] [19] and [20] We selected stars that have solutions 
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Fig. 18. Comparison between the temperature determined from 

photometry (our method) and spectroscopy (compiled catalog). 

Dots are all stars independent of the weight of the solution and 

the rms-error of the estimated temperature, crosses are stars with 

a weight of the solution w > 0.8, photometric errors better than 

0.05 mag in all bands, cr. ^pi, < 0.05 and Av[inag]/d[knc] e 

^ <■ / / , , 

[0.1,2.0] (criterion 1 & 3 as described in Sect. l4T2b . The right 

panel shows the corresponding histograms for both subsamples. 
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Fig. 19. Comparison between gravities determined from pho- 
tometry (our method) and spectroscopy (compiled catalog). Dots 
are all stars independent of the weight of the solution and the 
rms-error of the estimated gravity, crosses are stars with a weight 
of the solution w > 0.8, photometric errors better than 0.05 mag 
in all bands, (Xiog^,,/, < 0.5 and Av[mag]/d[kpc] e [0.1,2.0] (cri- 
terion 1 & 3 as described in Sect. 14.3.2b . The right panel shows 
the corresponding histograms for both subsamples. 



with weight w„ > 0.8, uncertainties of the input photometry bet- 
ter that 0.05 mag, rms-errors of the estimated parameters (disper- 
sion of the points around the weighted average parameter for the 
solution) better that 0.05 dex for log T^/f, 0.5 dex for log g, and 
0.3 dex for [M/H]. It becomes clear that gravity is the least ac- 
curate unknown in estimating of astrophysical parameters: in the 
case of calibration catalog we were able to estimate gravity only 
for 2 % of the stars with an appropriate precision (better than 0.5 
dex). We note that only 1186 stars from the calibration catalog 
have photometric errors better than 0.05 mag and only 21 stars 
have photometric errors better than 0.02 mag in all bands. We 
discuss the reasons for failures in the estimation of parameters 
below. We conclude that the combination of criteria 1 & 3 works 
best in selecting stars with the best estimation of astrophysical 
parameters. 

5.1. The problem of the selection of the synthetic model grid 

There are a number of libraries of synthetic models of atmo- 
spheres is available today. The most popular among them are 
the BaSeL library based on Kurucz models with an addition of 
synthetic and semi-empirical models for cool stars (Lejeune et 

1997 ), original Kurucz models (ATLAS9 program, Castelli 
& Kurucz. l20fel l. and MARCS models (Gustafsson et al.. l2003l l. 
Even the models of the same author for the same parameters (mi- 
croturbulent velocity etc.) differ due to improvements that the 
authors introduced into different versions of models. 

We do not discuss the different choices of theoretical models 
leaving this subject to the user of the method, but we note that 
some of the discrepancies in the parameters estimated via our 
method in comparison with the parameters from spectroscopic 
estimations can be due to the different model sets used for pho- 
tometric (performed by us) and spectroscopic (represented in the 
calibration catalog) estimations of astrophysical parameters. 

The direct comparison of Kurucz models with the SEDs of 
real stars (see, for example, Straizys et al., 1 19971 ) shows prob- 
lems with the ultraviolet part of spectra for main-sequence stars, 
as well as an inconsistency in the synthetic SEDs with real SEDs 



for late type stars (K7-M stars). This will not influence our re- 
sults for real stars as the lowest temperature estimation in Fig. [18] 
turned out to be 4200 K but must be taken into consideration in 
future work as late type stars are in the area with the best esti- 
mation of parameters possible (see Figs. [T3] and [T4b. 

Both photometric and spectroscopic determinations of astro- 
physical parameters based on synthetic models suffer from the 
uncertainties and systematic errors of the selected synthetic grid. 
We urge the users of this method to always make a reference to 
the synthetic models used and to accept the resulting parameters 
as parameters in the system of the selected synthetic models. A 
detailed study of the systematic errors of astrophysical param- 
eters due to problems of synthetic models is beyond the scope 
of this paper. There are many reasons for such systematics, for 
example, the convection treatment adopted for models and the 
influence of models on fundamental parameters is discussed in 
Heiter et al. ( 120021 ), who estimate that the systematic effect in 
temperature could be up to 400 K. 

The estimation of the systematics in the photometry 
due to the difference between synthetic models (Martins & 
Coelho, 12007! ) shows that, for Johnson B, V, J, H, and K bands, 
they are within 0.05 mag, which can be accepted as an upper 
limit for this type of uncertainties in photometry. The rough es- 
timation of the influence of these uncertainties on the resulting 
parameters can be done in the same way as the estimation of 
biases due to errors in the input photometry ( Figs. [T3] and [T4li. 

5.2. The problem of stars outside the grid 

As we can see from the previous section and Figs. [T5l[T6l[T7] we 
reached quite good agreement between spectroscopic and pho- 
tometric estimation of astrophysical parameters, however there 
are some cases of a large differences between our estimated pa- 
rameters and parameters from spectroscopy. To understand the 
reason for these differences we searched the Simbad database 
for the information about each star with "true" (spectroscopic) 
parameter out of the 3-cr range from the estimated parameters. 
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Fig. 20. Comparison between metallicities determined from pho- 
tometry (our method) and spectroscopy (compiled catalog). Dots 
are all stars independent of the weight of the solution and the 
rms-error of the estimated metallicity, crosses are stars with a 
weight of the solution w > 0.8, photometric errors better than 
0.05 mag in all bands, ct^m/h] < 0.3 and Av[mag]/d[kpc] e 
[0.1,2.0] (criterion 1 & 3 as described in Sect. l4T2b . The right 
panel shows the corresponding histograms for both subsamples. 



Table |4] summarizes the results and indicates the reason for the 
difference between them. 

The two stars with the greatest differences in the estimation 
of temperature are stars with peculiarities in the spectra. One of 
the cases where we failed to determine gravity is a Cepheid. For 
estimating the metallicity, we failed in the case of two close bina- 
ries, an RS CVn-like variable and a Cepheid (HD 101602, both 
gravity and metallicity were estimated wrong for this variable 
star). 

There is only one star for which the reason for the failure is 
unclear: BD -1-30 261 1. This star is a suspected Hipparcos binary 
star, although the radial velocity measurements did not prove the 
binarity (Sperauskas & Bartkevicius ,2002) . 

In all cases of failure, the stars in question did not satisfy the 
assumptions underlying our classification procedure. Binaries 
and variables cannot be estimated using the synthetic grid of 
models for "normal" single stars. Exotic stars (from the point 
of view of chemical abundance) like carbon stars were rejected 
at the very first step of the method by the method itself, which 
failed to find satisfying synthetic grid points. The same happens 
if we try to use photometric observations of a galaxy as an input 
(we used the catalog of Gavazzi & Boselli, [19961 ). In the case of 
binaries and variable stars we may get a wrong estimation of the 
parameters, if we have no idea about the true nature of the star 
and if we have no information about the distance of the star and 
the extinction. 



6. Results and conclusions 

We developed and tested a new method determining the astro- 
physical parameters from broad-band photometry. This method 
does not rely on a predefined model for the stellar population, 
but requires assuming the law for the extinction in the direction 
of the estimated star. In principle, the extinction law can be de- 
termined by the method itself with the use of stars with a unique 
solution for their astrophysical parameters, which is independent 
of the size of extinction. 



We used an interval-cluster analysis as the mathematical core 
for the method. For a given input photometric system (here 
Johnson B,V, and 2MASS J, H, and Ks), this method allows us 
to select the regions in the parameter space where a (unique) so- 
lution is possible and to identify regions in the parameter space 
where no solutions are possible due to the complicated topology 
of the color and parameter spaces. The astrophysical parameters 
determined with the use of this method are defined in the system 
of the selected models (synthetic or empirical) and suffer from 
all the systematic errors which are inherent in these models. 

We have compiled a catalog of stars with astrophysical pa- 
rameters known from spectroscopic observations to calibrate the 
Kurucz models and the estimated parameters. 

In the Johnson B,V, and 2MASS J, H, and Ks photomet- 
ric system and with mean errors of input colors better than 
0.01 mag the best results are achieved for main sequence G- 
K stars (most preferable region, G6 to K9). In this case the 
parameters are accurate to criogj-^^^ < 0.05, criogg < 0.3, and 
o"[M/H] < 0.3. If the accuracy of the color measurements drops 
to 0.05 mag, the most preferable region with the accuracy de- 
fined above (criog 7-__^.^. < 0.05, o-jogg < 0.3, and cr[M/H] < 0.3) 
shrinks to main sequence K6-K9 stars; meanwhile, the accuracy 
for the rest of the most preferable region drops dramatically, es- 
pecially for gravity. Nevertheless, even for 0.05 mag the accu- 
racy of the determined parameters for the rest of most preferable 
regions (main sequence G6 to K5) stays within criogT^ff < 0.1, 
(Tiogg < 2.0, and <T[m/h] < 0.4. The definition of "the most 
preferable regions" depends on the combination of photometric 
bands, desired precision for the result and the precision of the 
input photometry. 

We find that the biggest problem in determining astrophys- 
ical parameters is not to determine a solution but to distinguish 
between possible solutions that correspond to the input colors. 
We propose and tested a simple criterium for the selection of a 
solution based on the estimated extinction. Final tests with the 
compiled catalog proves the reliability and effectiveness of the 
method by comparing spectroscopic (literature) and photomet- 
ric (our method) determinations of parameters. 

The proposed method was initially designed to search for 
sparse stellar groups in deep sky surveys, but it also suits the 
study of populations in the Galaxy without a predefined model 



Table 4. Reason for the failure of the parameter estimation 



log T,ff 


name 




log Tl"^ 




HD 18078 


4.003 


4.654 


AOp 


HD 96446 


4.401 


4.635 


B2IIIp 




name 




log gP" 




BD +30 2611 


1.13 


3.20 




HD 101602 


1.40 


4.88 


UZ Cen, 








Cep-type variable 


[M/H] 


name 


[M/HY" 


[M/H]"" 




BD +30 2611 


-1.40 


-0.33 




HD 101602 


0.14 


-2.15 


UZ Cen, 








Cep-type variable 


HD 154338 


0.0 


-1.33 


V991 Sco, 








RS CVn-type variable 


HD 3940 


-0.39 


-1.85 


V755 Cas, 








Algol-type 


HD 6582 


-0.73 


-1.93 


SB 
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for the distribution of stars. It is possible to use the method to 
determine an extinction law. 
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